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ABSTRACT 

We present Herschel/HIFI observations of 30 transitions of water isotopologues toward the high- 
mass star forming region NGC 6334 I. The line profiles of H^O, Hi/O, H^ 8 0, and HDO show a 
complex pattern of emission and absorption components associated with the embedded hot cores, a 
lower-density envelope, two outflow components, and several foreground clouds, some associated with 
the NGC 6334 complex, others seen in projection against the strong continuum background of the 
source. Our analysis reveals an H2O ortho/para ratio of 3 ±0.5 in the foreground clouds, as well as the 
outflow. The water abundance varies from ~ 10~ 8 in the foreground clouds and the outer envelope to 
~ 10~ 6 in the hot core. The hot core abundance is two orders of magnitude below the chemical model 
predictions for dense, warm gas, but within the range of values found in other Herschel/HIFI studies 
of hot cores and hot corinos. This may be related to the relatively low gas and dust temperature 
(~ 100 K), or time dependent effects, resulting in a significant fraction of water molecules still locked 
up in dust grain mantles. The HDO/H2O ratio in NGC 6334 I, ~ 2 x 10~ 4 , is also relatively low, but 
within the range found in other high-mass star forming regions. 
Subject headings: stars: formation — ISM: molecules 
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1. INTRODUCTION 

Water is predicted to be one of the most abundant 
species in the interstellar medium (ISM) and, besides 
CO, the main reservoir of oxygen. In warm regions 
(Tdust > 100 K), where water sublimates from icy grain 
mantles, or is efficiently formed in the gas phase through 
reactions involving atomic oxygen, its gas-phase abun- 
dance may be as high 10 -4 with respect to H2 (e.g., Cec- 
carelli et al. 1996; Doty et al. 2002; Rodgers & Charnley 
2003), although in the presence of strong UV radiation 
the water abundance is typically much lower, ~ 10~ 7 , 
due to photodissociation (Hollenbach et al. 2009). Con- 
sequently, water may be an important coolant, espe- 
cially in high-mass star forming regions, where high gas 
and dust temperatures (> 100 K) and high densities 
(> 10 6 cm~ 3 ) are common. At lower temperatures 
(Tdust < 100 K), water molecules freeze-out onto dust 
grains and the resulting gas-phase H2O abundance is 
typically lower, by several orders of magnitude, than the 
chemical model predictions for cold clouds (e.g., Boon- 
man et al. 2003; Doty et al. 2002). The freeze-out of 
water in cold clouds has been confirmed by the detection 
of water ice (e.g., Willner et al. 1982; Keane et al. 2001; 
Gibb & Whittet 2002). Because of the potentially impor- 
tant role of water in the energy balance in star forming 
regions, and therefore its influence on the collapse and 
fragmentation of cloud cores, its abundance is an im- 
portant parameter in understanding the star formation 
process. However, because of the complex interplay be- 
tween different routes of the gas chemistry, grain-surface 
chemistry, freeze-out and desorption mechanisms, as well 



observational challenges, the chemistry of water in the 
interstellar medium (ISM) is not yet well understood. 

Most low-excitation H2O lines, especially the ground 
state ortho and para transitions, are not observable from 
the ground, because of the low atmospheric transmission 
due to the terrestrial water vapor. Space missions have 
already provided ample water observations in various 
ISM environments (van Dishoeck & Helmich 1996 — ISO; 
Melnick & Bergin 2005— SWAS; Bjerkeli et al. 2009— 
Odin). However, these have been limited to a small 
number of transitions (SWAS, Odin), low spectral res- 
olution (ISO), and low spatial resolution (> 1'). Pre- 
vious observations of high-mass star forming regions in- 
dicate an H2O abundance as high as ~ 10~ 4 in warm 
(Tdust > 100 K) regions, whereas in dense environments 
at temperatures below 100 K, where freeze-out takes 
place, the H2O abundance can be as low as 10 -8 (e.g., 
Boonman et al. 2003; Snell et al. 2000; Bergin & Snell 
2002; Herpin et al. 2012). The water abundance in low- 
mass star forming regions ranges from 5 x 10~ 7 in the 
envelopes to ~ 10~ 6 in the hot corinos (Ceccarelli et al. 
2000) and 10~ 5 in the outflows of Herbig Haro objects 
(Liseau et al. 1996). In diffuse clouds water abundances 
between 10 -7 and 10 -6 have been derived (Cernicharo 
et al. 1997; Moneti et al. 2001). 

Here we present water observations carried out with 
the HIFI instrument (de Graauw et al. 2010) aboard 
the Herschel Space Observatory (Pilbratt et al. 2010), 
which has provided for the first time the ability to ob- 
serve multiple water lines with high spectral and spatial 
resolution. Using this instrument, several investigations 
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of water in star forming regions have already been car- 
ried out. Mclnick et al. (2010) derived a water abun- 
dance of 1 — 7 x 10~ 5 in Orion KL, in agreement with 
chemical model predictions for warm gas. The Herschel 
guaranteed time key program "Water In Star- Forming re- 
gions with Herschel" (WISH; Van Dishoeck et al. 2011) 
is dedicated to investigations of water in some 80 galactic 
sources. Early WISH results show that the H2O abun- 
dance in the shocked gas in NGC 1333 is similar to the 
abundance in Orion KL (Kristensen et al. 2010), whereas 
Caselli et al. (2012) derived a much lower water abun- 
dance of > 10~ 9 within the central few thousand astro- 
nomical units in the prestellar core LI 544. Chavarria 
et al. (2010) investigated water in the high-mass proto- 
star W3 IRS5 and derived an H^O abundance of 10 -4 
in the hot core and 10~ 9 — 10~ 8 in the cooler envelope. 
All these measurements can be explained by theoretical 
models. However, van der Tak et al. (2010) derived an 
H2O abundance of only 2 x 10~ 10 in the dense core and 
4 x 10~ 9 in the cooler envelope of the high-mass star- 
forming region DR 21. Flagey et al. (2013) derived a 
water abundance of 5 x 10~ 8 with respect to H2 in the 
diffuse clouds in the Galactic disk, while Sonnentrucker 
et al. (2013) found a factor of 3 enhancement in the wa- 
ter abundance in the material within 200 pc from the 
Galactic center, as compared to the disk. These strong 
variations in the derived water abundance demonstrate 
the importance of fully understanding the water chem- 
istry in a wide variety of interstellar environments. 

The abundance of HDO and the D/H ratio in water 
is another key parameter, because it is an indicator of 
the formation temperature. From the cosmic D/H ra- 
tio (1.5 x 10 -5 ; Oliveira et al. 2003) one would expect 
an HDO/H2O ratio of 3 x 10~ 5 . However, strong en- 
hancements in abundances of deuterated molecules, due 
to their lower zero-point vibrational energies as compared 
to the normal variants (Phillips & Vastel 2003), take 
place at temperatures < 50 K (Roberts & Millar 2000). 
Since water is largely frozen-out at these low tempera- 
tures, ices on dust grains contain a high fraction (0.2%- 
2%) of deuterated water (Dartois et al. 2003; Parise 
et al. 2003). HDO has been observed in Sgr B2(M) 
and (N) (Comito et al. 2003, 2010) with abundances 
varying from 2.5 x lO" 11 {T gas < 100 K) to 3.5 x 10~ 9 
(T gas > 200 KjJ and HDO/H 2 0- 10~ 3 . In the low- 
mass protostar L1448-mm a very high HDO/H2O ratio 
of 0.04 was found (Codella et al. 2010), which agrees 
with the abundance ratios in other low-mass protostars 
such as IRAS 16293-2422 (0.2%-3%; Parise et al. 2005) 
and NGC 1333 IRAS2 (> 1%; Liu et al. 2011). Obser- 
vations of sites of high-mass star formation, which show 
in general lower abundances of deuterated molecules, re- 
vealed HDO/H2O ratios between 1 x 10" 4 (W3 IRS5; 
van der Tak et al. 2006) and 0.02 (Orion KL; Bergin et 
al. 2010). The high abundance of HDO in these sources 
indicates that water molecules have been predominately 
formed at a cold, earlier evolutionary stage and have just 
recently been evaporated from dust grains into the gas 
phase, where they can be studied by means of submil- 
limcter spectroscopy. 



In this paper, we analyze spectra of water isotopo- 
logues toward NGC 6334 I, a relatively nearby (1.7 kpc; 
Neckel 1978) high-mass star forming region, which har- 
bors sites of many stages of protostellar evolution (Straw 
& Hyland 1989). Single-dish submillimeter continuum 
observations indicate a total gas mass of 200 M Q asso- 
ciated with NGC 6334 I (Sandell 2000). The molecular 
hot core emission has been studied extensively over the 
past decade (e.g., Beuther et al. 2008, 2007; Hunter et al. 
2006). The source is associated with an ultra compact 
HII region (de Pree et al. 1995) and shows a very line-rich 
spectrum (Schilke et al. 2006; Thorwirth et al. 2003). 
Furthermore, a bipolar outflow (Leurini et al. 2006; 
Beuther et al. 2008) and H 2 0, OH, CH 3 OH class II, and 
NH3 masers have been detected (e.g., Kraemer & Jack- 
son 1995; Ellingsen et al. 1996; Walsh et al. 2007). Rolffs 
et al. (2011) studied the radial density structure of sev- 
eral high-mass star forming cores, including NGC 6334 I, 
mainly based on single-dish HCN observations with the 
APEX telescope, and approximated the source structure 
as a centrally heated sphere with a power-law density 
distribution. In the case of NGC 6334 I, this leads to a 
density law n ~ r -15 , which describes well the envelope 
traced by the single-dish observations. However, inter- 
ferometric data (SMA) reveal that the internal structure 
of NGC 6334 I is much more complex — the hot core con- 
sists of four compact condensations within a 10" region 
that contain about 50% of the single-dish continuum flux 
(Hunter et al. 2006). High-resolution molecular observa- 
tions, combined with 3-d radiative transfer models, are 
required to fully understand the distribution of water in 
NGC 6334 I, however, simple models like the ones pre- 
sented here already provide a good starting point for a 
more sophisticated analysis of the HIFI data. 

The data presented here are taken from an unbiased 
HIFI spectral line survey, carried out as part of the Her- 
schel guaranteed time key program "Chemical HErschel 
Surveys of Star forming regions" (CHESS; Ceccarelli et 
al. 2010). In an earlier paper (Emprechtinger et al. 2010, 
hereafter paper I), we presented a preliminary analy- 
sis of twelve water lines in NGC 6334 I, including the 
In — 0oo and lio — loi fundamental rotational transi- 
tions, observed as a part of the Herschel science demon- 
stration program (SDP). Assuming that all ortho and 
para molecules are in their respective ground states, we 
determined an ortho/para ratio of ~ 1.6 ± 1 (3er uncer- 
tainty) in the foreground clouds, which appeared lower 
than the statistical value of 3, as seen for example in the 
outflow. The water abundance ranges between 10 -8 and 
4 x 10~ 5 , in agreement with theoretical models and pre- 
vious findings (e.g. Rodgers & Charnley 2003; Doty et 
al. 2002; Ceccarelli et al. 1996; Herpin et al. 2012)'. In 
addition, the H^O/H^O ratio of 3.7 ± 0.6 agrees well 
with the isotopic 18 0/ 17 ratio derived from formalde- 
hyde measurements (Wilson & Rood 1994). 

Here we present the complete HIFI data set of water 
observations toward NGC 6334 I, including 30 lines of 
H^ 6 0, H^ 8 0, H^O, and HDO. In the following sections, 
we present the observations (§2), along with a simple 
LTE analysis (§3), and full radiative transfer modeling 
of H;j, 6 lines (§ 4). We summarize the results in §5. 



1 At the molecular hydrogen densities typical of dense interstellar 
clouds, the gas and dust are usually well coupled, T gas = T^ ust . 



Water in NGC 6334 I 



3 



2. OBSERVATIONS 

HIFI observations of NGC 6334 I (a J20 oo = 
17 h 20 m 53.32 s , S J200Q = 35°46'58.5") presented here were 
carried out between 2010 February 28 and October 14, 
in the double beam switch mode (180" chopper throw). 
The data have been reduced using the Herschel Inter- 
active Processing Environment (HIPE; Ott 2010) ver- 
sion 5.1. The spectral resolution of the double sideband 
(DSB) spectra, obtained using the HIFI WBS spectrom- 
eter, is ~ 1.1 MHz. The DSB spectra have been observed 
with a redundancy of 8 for the SIS mixer bands (8 LO 
settings per IF bandwidth) and 4 for the HEB bands, 
which allows for the deconvolution and isolation of the 
single sideband (SSB) spectra (Comito & Schilke 2002). 
The deconvolved SSB spectra have been exported to the 
FITS format for subsequent analysis using the IRAM 
GILDAS packagfl Spectra presented here are equally 
weighted averages of the H and V polarizations, corrected 
for the corresponding aperture efficiencies (Roelfsema et 
al. 2012; Table 1). 

In the CHESS spectral survey, we detect sixteen H 2 6 0, 
six H^O, four H^O, and four HDO lines. Several other 
lines, especially of the rarer isotopologues, are blended 
with emission lines of other molecules, and thus provide 
only upper limits for the analysis. In Table 1 we sum- 
marize the observing parameters of the detected water 
lines. 

3. RESULTS & ANALYSIS 
3.1. Line Shapes 

Spectra of all 30 water transitions observed in 
NGC 6334 I are displayed in Figure 1. Their line shapes 
differ significantly, indicating that different lines trace 
different regions within the beam. From the line shapes, 
we can identify four major components. The first com- 
ponent, seen in many transitions, has a Gaussian line 
shape at a velocity of approximately -6.4 kms -1 , with a 
line width of approximately 5 kms -1 . This feature cor- 
responds to the envelope of NGC 6334 I. In the ground 
state transitions of H^O, but also in the H^O 2 2 i — 2i 2 
line, this component is seen in absorption against the 
background continuum; in most other transitions it is 
seen in emission. The second component is identified 
with one of the embedded cores. Some emission lines, 
especially those of highly excited lines of H 2 6 and lines 
of rare isotopologues, have central velocities close to - 
8 kms -1 , corresponding to the local standard of rest 
(LSR) velocity of the embedded core SMA 2 (Beuther 
et al. 2007). Thus we associate these lines with the 
dense, compact core. 

The third major component are the foreground clouds, 
some of which reside in the immediate vicinity of 
NGC 6334 I, while others are not physically associated 
with the NGC 6334 complex, van der Wiel et al. (2010) 
identified four foreground components at -3, 0, +6.5, and 
+8 kms -1 based on HIFI observations of CH. Three of 
these components (0, +6.5, and +8 kms -1 ) are also de- 
tected in the ground state H 2 6 transitions. Whereas 
the kms -1 component most likely originates in the lo- 
cal environment of NGC 6334 I, the features at +6.5 and 

2 http: / /www.iram.fr /IRAMFR / GILDAS / 



+8 kms -1 are associated with clouds in the Sagittarius 
spiral arm. 

The fourth major spectral component are the outflows, 
evidence of which is seen in almost all H 2 6 lines with 
J u < 3, as well as all ground-state transitions of the 
rarer isotopologues. In some transitions the outflows are 
detected in emission, in others in absorption, and in some 
cases P-Cygni profiles are seen. 

The physical properties of these four physical compo- 
nents are summarized in Table 2. In the following sec- 
tions we investigate in detail three of these components: 
the dense core, the outflow, and the foreground clouds, 
individually. Table 3 lists all transitions observed and 
specifies which components they trace. 

3.2. Foreground Clouds 
3.2.1. Ortho/Para Ratio Revisited 

Foreground absorption features are best seen in the 
ground state transitions of para-H 2 6 (In — 0oo) and 
ortho-H 2 6 (lio — lot, %12 — loi ), which are displayed 
in Figure 2. In paper I, we determined the ortho/para 
ratio in the foreground components using the lio — loi 
and the In — 0oo lines. Under the assumptions that all 
ortho and para water molecules are in their respective 
ground states, these observations yield a low ortho/para 
ratio of 1.6 + 1 (3er uncertainty), below the statistical 
value of 3. With the additional observations of the sec- 
ond ground state ortho-water line, the 2i 2 — loi transi- 
tion at 1669 GHz, we can independently constrain the 
excitation temperature, assuming a level population ac- 
cording to the Boltzmann distribution and determine the 
ortho/para ratio more reliably. In order to derive the 
ortho/para ratio, we first resample the spectra to a com- 
mon velocity resolution of 0.27 kms -1 , which allows us 
to compare column densities in individual velocity bins. 
The optical depth of the ground state transitions is cal- 
culated as 

lr ( T B- MTex) \ (1) 

where Tb is the measured brightness temperature, 
Tback is the corresponding background temperature, and 
3 v {T ex ) = h^ji^vlkT^ _y )S& the Ray i e igh- Jeans equiv- 
alent radiation field due to the molecular emission. We 
neglect the influence of the cosmic microwave background 
(CMB), which is negligible at the frequencies considered 
here (< 1.5 mK). In the case of the In — 0oo and 2i 2 — loi 
lines, the background radiation is well defined by the 
continuum level at nearby frequencies. For the lio — loi 
line, however, the background is dominated by the high- 
velocity emission from the outflow. We thus applied a 
linear baseline fit in the vicinity of the absorption com- 
ponents in order to determine the correct background 
level for the optical depth computations!! Optical depths 
less than 2.5 are determined quite reliably, however, for 
t > 2.5, the uncertainty of r is comparable to r. Thus 
we do not include these data points in the subsequent 

3 We assume that the absorption completely covers the contin- 
uum. This is a reasonable assumption, as for example the OH data 
of Brooks & Whiteoak (2001) show that the 6 kms" 1 absorption is 
seen not only toward NGC 6334, but also toward other nearby HII 
regions. In addition, the deepest absorption seem in the ground 
state water lines is very close to zero. 
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I — to *~ ^ to 

Fig. 1.— HIFI spectra of the 30 transitions of H^ e O, H^ 7 0, H^ s O, and HDO toward NGC 6334 I. The intensity scale is brightness 
temperature, in K, corrected for the aperture efficiency. Dashed vertical lines mark the velocity of the dense core SMA 2 at —8 kms -1 . 
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TABLE 1 

Water transitions observed in NGC 6334 I. 



I sotopologue 


Transition 


HIFI Band 


Frequency 


E up 


Beam Size 


Vap 










(GHz) 


(K) 


(") 




p-H«0 


6'>4 




1A 


488.491 


867.6 


43.4 


0.68 




2i i 




2B 


752.033 


137.0 


28.2 


0.67 




4oo 
^-ZZ 




3B 


916.172 


454.5 


23.2 


0.67 






- 4 3 i * 


4A 


970.315 


599.1 


21.9 


0.65 




2o2 


- In 


4A 


987.927 


100.9 


21.5 


0.65 




111 


- Ooo * 


4B 


1113.343 


53.5 


19.0 


0.65 




422 


— 4l3 


5A 


1207.639 


454.5 


17.6 


0.55 




220 


- 2n 


5A 


1228.789 


196.0 


17.3 


0.55 


o-H^O 


lio 


- I01 * 


1A 


556.936 


61.0 


38.1 


0.68 




532 


- 4 4 i 


IB 


620.701 


732.4 


34.2 


0.67 




3i2 


— 3o3 


4B 


1097.365 


249.5 


19.3 


0.65 




3l2 


- 2 2 1 * 


5A 


1153.127 


249.5 


18.4 


0.55 




634 


— 541 


5A 


1158.323 


934.2 


18.3 


0.55 




321 


— 3i2 


5A 


1162.912 


305.4 


18.2 


0.55 




221 


-2l2 


6B 


1661.008 


194.2 


12.8 


0.63 




2l2 


-loi 


6B 


1669.905 


114.4 


12.7 


0.63 


p-H^O 


2n 


— 2o2 


2B 


745.320 


136.4 


28.4 


0.67 




202 


- Ill 


4A 


994.675 


100.7 


21.3 


0.65 




111 


- Ooo * 


4B 


1101.698 


52.9 


19.3 


0.65 


o-H^O 


lio 


-loi 


1A 


547.676 


60.5 


38.8 


0.68 




3l2 


— 303 


4B 


1095.627 


248.9 


19.4 


0.65 




2l2 


- loi 


6B 


1655.867 


113.7 


12.8 


0.63 


p-H^'O 


2n 


— 2 ( )2 


2B 


748.458 


136.7 


28.3 


0.67 




In 


- Ooo * 


4B 


1107.167 


53.2 


19.2 


0.65 


o-H^O 


lio 


-loi * 


1A 


552.020 


60.7 


38.4 


0.68 




3l2 


- 3() 3 


4B 


1096.414 


249.2 


19.3 


0.65 


HDO 


2n 


— 2()2 


IB 


599.927 


95.3 


35.3 


0.67 




In 


- Ooo 


3B 


893.639 


43.9 


23.7 


0.67 




2o2 


- 1(11 


3B 


919.311 


66.5 


23.1 


0.67 




2n 


- lio 


4A 


1009.945 


95.3 


21.0 


0.65 



Notes — Entries in the tabic are: the isotopologue and transition, the HIFI mixer band, the rest frequency (GHz), the upper state energy 
(K), the Herschel beam size (arcsec), and the aperture efficiency (Roelfsema et al. 2012). An asterisk marks lines that have already been 
presented in paper I. 



TABLE 2 

Physical properties of the gas components in NGC 6334 I. 





Envelope 


Hot Core 


Outflow 




Foreground 










Broad 


Narrow 


1 


2 3 


4 


fLSR (kms 1 ) 


-6.4 


-8.2 


-10 


-10 


-3 


+6.5 


+8 


Av FWHM (kms" 1 ) 


5.4 


4.0 


80 


32 


5.0 


3.0 1.2 


2.5 


T gas (K) 


20-100 


> 100 


> 50 


130 




~ 30 




n(H 2 ) (cm- 3 ) 


5 x 10 4 -5 x 10 6 


10 7 -10 8 




2 x 10 7 




< 10 4 
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TABLE 3 

Spectral components seen in the water spectra. 



Isotopologue 


Transition 


Frequency 


Main 


Dense 


Outflow 


F oreground 








(GHz) 


C OTTipOTlCTlt 








D-Hi e O 


624 


- 7 17 


448.5 


N 


Y 


N 


N 




2n 


— 2o2 


752.0 


E 


N 


Y 


n; 




422 


- 331 


916.2 


N 


Y 


N 


N 




524 


- 4 3 i 


970.3 


N 


Y 


N 


N 




202 


- In 


987.9 


E 


N 


Y 


n; 




111 


- Ooo 


1113.3 


A 


N 


pi> 


Y 




422 


- 4l3 


1207.6 


N 


Y 


N 


N 




220 


- 2n 


1228.8 


E 


N 


p 


N 


o-H^O 


lio 


- I01 


556.9 


A 


N 


yb 


Y 




532 


- 44i 


620.7 


N 


Y 


N 


N 




3l2 


— 3o3 


1097.4 


E 


N 


p 


N 




3l2 


— 221 


1153.1 


E 


N 


Y 


K 




634 


— 541 


1158.1 


N 


Y 


N 


N 




021 


— 912 


1162.9 


E 


N 


N 


N 




9^ 
^21 


— ^12 


1661.0 


A 


N 


P 


N 




9i ^ 
^12 


1 

— J-01 


1669.9 


A 


N 


p6 


Y 


p-H^O 




2ll 




— 202 


745.3 


N 


Y 


N 


N 




^02 


— J-ll 


994.7 


N 


Y 


N 


N 




J-ll 


— Uoo 


1101.7 


N 


N 


A b 


N 


o-H^O 


410 


— 101 


547.7 


N 


N 


A b 


N 




3l2 


— 303 


1095.6 


N 


Y 


N 


N 




2l2 


-loi 


1655.9 


N 


N 


A» 


N 


p-H^O 


2n 


— 2o2 


748.5 


N 


Y 


N 


N 




111 


~Ooo 


1107.2 


N 


N 


A b 


N 


o-H^O 


lio 


-loi 


552.0 


N 


N 


A b 


N 




3l2 


- 3o3 


1096.4 


N 


Y 


N 


N 


HDO 


2ii 


— 2q2 


599.9 


N 


Y 


N 


N 




In 


- Ooo 


893.6 


N 


N 


A b 


N 




2o2 


-loi 


919.3 


N 


Y 


N 


N 




2n 


- lio 


1009.9 


N 


Y 


N 


N 



Notes — Main Component: E= seen in emission, A=seen in absorption. Dense Core: Y=a component detected at a velocity of —8 kms -1 , 
and thus associated with the dense cores. Outflow: Y=outflow detected in emission, A=bluc lobe of the outflow in absorption, P= P-Cygni 
profile. The superscript b indicates that the broad outflow component is seen. Foreground: Y=foreground components are detected. N in 
all three columns means that this component does not show up in the respective transition. 



analysis. From the optical depth, we determine the pop- 
ulation in the ground state per velocity interval as 

_ 8ir gi t 

'ground- x s Agul _ e _ hv/kTex W 

where A is the wavelength, A the Einstein spontaneous 
emission coefficient, gi and g u the statistical weights for 
the lower and the upper state, respectively, and T ex the 
excitation temperature. 

Combining equations (1) and (2), we obtain for each 
transition an equation with two unknowns, N grounc i and 
T ex . Assuming a level population according to the Boltz- 
mann distribution implies the same excitation tempera- 
ture for all three transitions detected in the foreground 
clouds. Furthermore, the population of the ground state 
for the two ortho-H 2 lines is identical, as they both 
connect to the loi level. We can thus eliminate it from 
the equations and compute the excitation temperature 
directly. This leads to T ex = 6.5 K. 

As a consistency check, we derive an upper limit for 
the population in the In para-H 2 state, and thus for 
the excitation temperature, from the non-detection of the 
foreground clouds in the 2o 2 — In transition. The upper 
limit on the optical depth in this line is 0.009 leading to 
an upper limit of T ex < 8.7 K, in agreement with the 
value derived from the ortho-H 2 lines. 



Equations (1) and (2), with the derived excitation tem- 
perature of 6.5 K, lead to an ortho/para ratio of 3±0.5 in 
the foreground clouds, in good agreement with the statis- 
tical value of 3. Figure 3 shows the calculated ortho/para 
ratio as a function of velocity. Comparing these results 
with the simple analysis in paper I (assuming that all 
population is in the respective ground state, leading to 
an ortho/para ratio of 1.6) reveals that the influence of 
the excitation temperature on the derived column den- 
sity of the para-H20 is negligible, due to the high energy 
of the first excited state (In, i?=53.4 K). Even the parti- 
tion function of ortho-H 2 changes by less than 3%. The 
main reason for the discrepancy between the two meth- 
ods is that the quantity J u (T ex ) at the frequency of the 
lio — loi ortho-H 2 line, 557 GHz, is a sensitive function 
of T ex . Thus the optical depth and therefore the column 
density of ortho-H 2 was underestimated in paper I by 
approximately a factor of 1.9, leading to the low value of 
1.6 for the ortho/para ratio. 

The assumption of a thermalized level population may 
not be fully justified, considering the high critical densi- 
ties of 6.4 x 10 7 cm- 3 and 1.12 x 10 8 cm" 3 of the l w - l i 
and 2i2 — loi transitions, respectively. However, since the 
critical density of the 2 i2 — loi line is higher, it is the ex- 
citation temperature of this transition, which would be 
overestimated by assuming LTE. But, as shown above, 
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Fig. 2. — Spectra of the three ground state transitions of ortho- 
and para-Hj O toward NGC 6334 I. The intensity scale is bright- 
ness temperature, in K, corrected for the aperture efficiency. 



only the lio — loi line is sensitive to changes in T ex . 
Therefore a deviation from the Boltzmann distribution 
does not have significant influence on the results pre- 
sented here. 

An excitation temperature of 6.5 K for the H2O is rel- 
atively high for diffuse clouds. However, some of the 
foreground components toward NGC 6334 1 are also seen 
in absorption in the ground state para-NH3 transitions, 
which indicates the presence of gas denser (~ 10 5 cm~ 3 ) 
than typically found in diffuse clouds. This analysis 
nicely illustrates the importance of the HIFI capabili- 
ties to observe multiple transitions for the determination 
of accurate column densities and molecular abundances. 

3.2.2. Water Abundance in the Foreground Clouds 

To determine the water abundance in the foreground 
clouds, we use the HF (1-0) absorption spectrum (Em- 
prcchtinger et al. 2012) to estimate the H2 column den- 
sity. HF has been established as an excellent tracer of H2 
(Neufeld et al. 2005; Neufeld et al. 2010; Sonnentrucker 
et al. 2010; Monje et al. 2011). Following these stud- 
ies, we use an HF abundance of 1 x 10~ 8 with respect to 
H2, as determined on several lines of sight. Because of 
the high Einstein A coefficient (2.423 x 10 2 s~ x ) and the 
low collisional excitation rates (Guillon et al. 2011), the 
critical density of HF (1-0) is ~ 6 x 10 9 cm~ 3 , and thus 
approximately three orders of magnitude higher than the 
critical density of the ortho-H20 lio — loi transition. 
Therefore, we assume that all HF molecules are in the 
ground state (see Emprechtinger et al. 2012). Unfortu- 
nately, the HF (1-0) line, like H^O , quickly becomes 
optically thick (r > 2.5) at the line center and therefore 
we can determine the HF column density only in the line 
wings, where we derive an average H2O abundance of 
1 — 3 x 10~ 8 in the foreground clouds (see Fig. 4), com- 
parable to typical values determined in low-temperature 
regions (T dust < 100 K). 




Fig. 3. — Water ortho/para ratio toward NGC 6334 I in the ve- 
locity range from -4 to 10 kms - , corresponding to the foreground 
clouds. Green line shows the P-H2O absorption (middle panel in 
Fig. 2). Red line indicates the statistical limit of 3. 



3.3. Outflows 

In paper I, we have already shown that NGC 6334 I 
harbors two outflow components. One of these compo- 
nents shows broad lines (from -85 to +65 kms -1 ) sim- 
ilar to the CO outflow detected by Leurini et al. (2006), 
who derived a gas temperature > 50 K and a density 
> 1 x 10 4 cm~ 3 at the peak position of the red lobe, and 
T gas > 15 K and n > 1 x 10 3 cm' 3 toward the blue lobe. 
This outflow component was already analyzed in paper 
I, where we determined that most of the water molecules 
are in the ground state, as this outflow component is not 
seen in any excited lines. This is in contrast to HIFI ob- 
servations of other high-mass star forming regions that 
do show bright high-velocity emission in excited water 
lines (e.g. W3 IRS5; Chavarn'a et al. 2010) and may 
be related to the relatively low density in the outflowing 
gas (Leurini et al. 2006), resulting in subthermal excita- 
tion of water molecules. Based on the H^O and H3, 8 
observations, we derived and optical depth of ~ 300 for 
the H2O In — 0oo transition in the velocity range -35 
to -20 kms -1 . A comparison with LVG models (Radex; 
van der Tak et al. 2007) shows that our findings are in 
agreement with the lower limits for the temperature and 
density found by Leurini et al. (2006). However, despite 
its high optical depth, this outflow component does not 
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TABLE 4 

Integrated line intensities of H3, 6 transitions detected 

IN THE NARROW OUTFLOW. 



x(H 2 0) 
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Fig. 4. — H2 column density (upper panel, based on HF absorp- 
tion data) and H2O abundance (lower panel) in the foreground 
clouds on the line of sight toward NGC 6334 I, plotted as a func- 
tion of the LSR velocity. 



absorb all background radiation, leading to the conclu- 
sion that only approximately 17% of the total continuum 
of NGC 6334 I is subject to absorption by the outflow 
(area coverage factor of 17%). The water abundance in 
the broad line outflow is 4 x 10" 5 (paper I), due to sub- 
limation or sputtering of water ices. The second outflow 
component is much narrower, ranging from -45 to +25 
kins" 1 . Furthermore, the narrow outflow is also seen in 
some excited lines, which indicates that this component 
is warmer and denser. It is not clear, if these two outflow 
components represent two separate outflows or different 
regions within one outflow. Within the measurement un- 
certainties, the central velocity of both outflow compo- 
nents is -10 kms" 1 , which is consistent with the latter 
scenario. This is further supported by interfcrometric 
observations of Beuther et al. (2008) who found only 
one outflow in their HCN (1-0) data. 

3.3.1. Warm and Dense Outflow 

The three ground state water transitions are not very 
well suited to investigate the warm and dense outflow 
component, because the much more massive, but less ex- 
cited broad outflow dominates the spectra. Furthermore, 
foreground absorption features hamper a clear view of 
the high-density outflow. However, the dense outflow 
features can be seen clearly in the moderately excited 
2o2 — In, 2n — 2 2, and 2 2 o — 2n para-H 2 lines, as well 
as the 3i2 — 3o3 ortho-H 2 line (Figure 5). 

In the 2 2 o — 2n and 3i 2 — 3o3 transitions, the warm 
and dense outflow shows a P-Cygni profile, but the blue 
wing is also weaker than the red wing in the remaining 
two transitions. The intensity difference between the red 
and blue wings is approximately 30% of the continuum 
in all four transitions, leading to the conclusion that all 
four lines are optically thick (r > 3). The lack of any 
corresponding outflow features in the H 2 8 lines gives 
an upper limit of 80 for the optical depth. 



Transition 


Frequency 


jT B dv 




(GHz) 


(K kms" 1 ) 


2o2 — 111 


987.9 


24.1 ±2.4 


2n — 2 ( )2 


752.0 


26.2 ±2.6 


220 — 2n 


1228.8 


10.2 ± 1.0 


3l2 — 3()3 


1097.4 


22.7 ±2.3 



The red lobe of the outflow, which is expected to re- 
side behind the continuum source (the blue wing is seen 
in absorption, while the red one is not) can be used to 
determine the physical conditions in the outflowing gas. 
The integrated intensities emitted by the warm and dense 
outflow are listed in Table 4. 

Using Radex (van der Tak et al. 2007), which cal- 
culates non-LTE radiative transfer, with the collisional 
rates from Faure et al. (2007), and assuming an isother- 
mal, homogeneous medium, we derive the temperature 
and density in the outflow. Input parameters for the 
Radex models are: the background radiation tempera- 
ture, the gas kinetic temperature, the H 2 density, the 
column density of water molecules, and the line width. 
For these parameters, the excitation temperature, op- 
tical depth, and radiation temperature of each line are 
calculated. 

We set the background temperature to 2.73 K and the 
line width to 30 kms" 1 , a value which we derive from a 
fit of several outflow features. Since we do not know the 
size, and thus the beam dilution of the outflow compo- 
nent, the observed line intensities do not provide good 
constraints. For this analysis we use an H 2 column 
density of 1 x 10 18 cm" 2 , resulting in an optical depth of 
~ 25, well within the limits 3 < t < 80. A x 2 analysis 
of the model calculations reveals that the observations 
are best reproduced by an H 2 density of 2.5 x 10 7 cm" 3 
and a kinetic temperature of 130 K. These results are in 
agreement with the non-detection of the dense outflow in 
the highly-excited lines. Whereas the H 2 density of the 
dense outflow component is quite well defined, with pos- 
sible solutions ranging from 1.5 x 10 7 to 5.0 x 10 7 cm" 3 , 
the temperature is less constrained and acceptable fits 
can be found for values between 60 K and 200 K. 

The water column density in the dense outflow is not 
well constrained either, since all four H^O lines are op- 
tically thick. Taking the limits for r given above and 
the best fit results for the density and temperature, we 
derive a water column density between 5.0 x 10 16 and 
1.2 x 10 18 cm" 2 . Using the upper limit for the H 2 col- 
umn density does not change the results for temperature 
and the H 2 density significantly. The best result for the 
lower limit on the H 2 column density can be found at a 
temperature of 94 K and an H 2 density of 3.6 x 10 8 cm" 2 . 
Taking the uncertainty due to the unknown H 2 column 
density into account, we derive a temperature between 
60 K and 200 K and an H 2 density higher than 2 x 10 7 
cm" 2 for the dense outflow component. 

3.3.2. HDO in the broad outflow 

Besides H^O, H^O, and H^O, which have been pre- 
sented in paper I, the HDO In — 0oo ground state tran- 
sition is detected in absorption as well. In Figure 6 we 



Water in NGC 6334 I 



9 



T B (K) 




-50 50 

Velocity (km/s) 




-50 50 

Velocity (km/s) 

Fig. 5. — Excited H2O lines showing the narrow outflow component 
the aperture efficiency. 

show a comparison of the HDO and H3> 7 line profiles. 
The absorption feature covers the velocity range from —6 
kms -1 to approximately -50 kms -1 and thus originates 
in the blue lobe of the broad outflow. 

To calculate the column density of HDO in the outflow, 
we use equations (1) and (2), assuming that all molecules 
are in the ground state. This assumption is justified by 
the non-detection of emission in the red lobe of the out- 
flow, which sets an upper limit for the HDO excitation 
temperature. 

The main difference between the calculation of the col- 
umn density here and in foreground clouds (§3.2) is that 
the broad outflow does not cover the entire continuum 
source (see paper I). From the absorption features of 
H^ 6 and H^ 8 f 10 - 00 , we determined that the H^ 6 
lio — 0oo transition is optically thick, but the absorp- 
tion feature lowers the intensity by only 17%. Thus we 
concluded in paper I, that only 17% of the continuum is 
covered by the outflow. To account for this low contin- 
uum coverage we have subtracted 83% of the continuum 
from Tb and Tb ac k in equation (1). 

This calculation leads to an HDO column density of 
1.17 ± 0.4 x 10 13 cm" 2 and a para-H^O column density 
of 1.21±0.4 x 10 13 cm" 2 . The HDO/p-Hf O ratio is thus 
1.0 ±0.5. Assuming an ortho/para ratio of three (see pa- 
per I), and an H^O/H^O ratio of 1600, the resulting 
HDO/H^ 6 ratio is 2.1 ±1 X 10~ 4 . This value is approx- 
imately an order of magnitude above the cosmic D/H 
ratio, but much lower than the HDO/H2O ratios found 
in low-mass protostellar cores (~ 10~ 2 ; e.g., Codella et 
al. 2010). 

A source of significant uncertainty in the HDO/H2O 
ratio is the assumption that all molecules are in their 
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respective ground states. An higher excitation tempera- 
ture would change the resulting HDO and H^O column 
densities. However, the influence of the excitation tem- 
perature on the abundance ratio is limited. The upper 
state energies of the two transitions differ by only 25%, 
and thus the relative population of the corresponding 
upper levels is very similar for any value of T ex . Fur- 
thermore the moderate optical depth of both lines limits 
the uncertainty due to J(T ex ) in equation (1). Overall, 
we estimate that the uncertainty of the HDO/H2O ratio, 
caused by the uncertainty in T ex is less than 20%. 

3.4. Dense Cores 

Hot, dense cores in NGC 6334 I have been observed 
interferometrically (Hunter et al. 2006; Beuther et al. 
2007, 2008) and their LSR velocities range from -7.6 to 
-8.1 kms -1 , which is 1.4 to 1.9 kms -1 lower than the 
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systemic velocity of the surrounding envelope. Therefore 
we use the velocity of the emission lines to identify the 
transitions originating from these dense cores. 

3.4.1. Hf O lines 

Five clearly detected and one tentatively detected 
H^O lines are identified to stem from the dense cores. 
Line parameters of these transitions are listed in Ta- 
ble 5. All five lines have upper state energies of 450 K or 
higher, indicating the presence of hot gas. Furthermore, 
none of these lines show any indication of an outflow, 
or a second component, suggesting that only the dense 
cores contribute to the observed emission. Unfortunately 
all but two of these lines are known to show maser ac- 
tivity in certain interstellar environments (Maercker et 
al. 2008; Neufcld & Melnick 1991). Therefore the lo- 
cal density and temperature structure, as well as the 
local radiation field, influence the observed line intensi- 
ties. None of these transitions show particularly narrow 
lines, and thus no compelling sign of maser activity, yet 
no solution assuming LTE can be found. The fact that 
the 624 — 7i7 H2O line, the line with the highest upper 
state energy, is much stronger than predicted by any LTE 
model with temperatures applicable to dense molecular 
clouds (< 1000 K) indicates that the excitation is clearly 
out of thermal equilibrium. We also attempted to model 
the line emission using Radex (§3.3.1), and with reason- 
able assumptions (kinetic temperature < 1000 K, n(H 2 ) 
from 10 5 to 10 s cm~ 3 ) could not find a satisfactory solu- 
tion for the hot core. The Radex models under-predict, 
similarly to the LTE models, the highest energy 624 — 7n 
transition by 1-2 orders of magnitude. Furthermore, for 
all lines except the 4 2 2 — 4i 3 transition, the calculated 
level population is either inverted, or the upper level is 
more populated than expected for reasonable values of 
the gas kinetic temperature. All these findings suggest 
that a full radiative transfer model, including pumping 
by the dust continuum radiation is necessary to fully ex- 
plain the H-^O emission from the hot core in NGC 6334 I 
(see §4). 

3.4.2. Rarer Isotopologues 

The dense core is not seen in any ground state transi- 
tion, because of the high optical depth of the surround- 
ing envelope material. However, moderately excited lines 
(J u = 2 or 3) of the rarer isotopologues H^O and H^O 
reveal the dense core emission. Contrary to the corre- 
sponding H^O lines, these moderately excited transi- 
tions of the rarer isotopologues do not show any sign of 
outflows or multiple components and the central veloc- 
ity of these lines confirms that they stem from the dense 
core (see Table 6). 

Figure 7 displays moderately-excited lines of H^O, 
H-^O, and H^O. Whereas the H^O line is approxi- 
mately a factor of ten stronger than the rarer isotopic 
lines, the corresponding H^O and H3; 7 lines have com- 
parable line strengths. The H^O/H^O ratios arc 1.21 
and 1.66 for the 2n — 2 02 line and the 3i2 — 3o3 line, re- 
spectively. This may be explained by an optically thick 
dense core embedded in a larger, optically thin envelope, 
only seen in H^O. Assuming similar excitation condi- 
tions for HJ^O and H^O, which is justified by the high 
densities expected in the dense core, and an H^O/H^O 



abundance ratio of 3.2 ± 0.2 (Wilson & Rood 1994) leads 
to t(H4 7 0) = 1.78 and 0.81 for the 2n-2 02 and 3i 2 -3 03 
lines, respectively. In addition to H^O and H^O, lines 
of HDO have also been detected (see Table 6). The cen- 
troid velocity of the HDO lines (-8.2 kms -1 ) confirms 
that these lines originate from the dense core, as well. 

We analyze the observed lines of all three rare water 
isotopologues, H^O, H^O, and HDO independently, us- 
ing XCLASS (Schilke et al. 1999; Comito et al. 2005), 
which calculates the level population assuming LTE, and 
from there the optical depth, J{T ex ), and Tg. The input 
parameters are the source size, excitation temperature, 
column density, velocity, and line width. LTE is in gen- 
eral not a valid assumption for water lines, due to their 
high critical densities (~ 10 7 cm -3 ), and thus any LTE 
analysis of water lines has to be treated with caution. In 
case of a non-LTE level population, both the temperature 
and the column density derived may be largely underesti- 
mated. However, as shown later by full radiative transfer 
modeling (see § 4) under the physical conditions present 
in NGC 6334 I, the water abundance is reproduced quite 
well assuming LTE. 

For our analysis we use a velocity of -8.2 kms -1 and 
a line width of 4 kms -1 . The resulting size of the dense 
cores is 4 ± 1", in good agreement with interfcrometric 
observations (Beuther et al. 2007; Hunter et al. 2006). 
The derived excitation temperatures and column densi- 
ties are listed in Table 7: the best fit for H^ 7 and HDO 
is 50 K, while a slightly lower value of 43 K is derived 
for H^O. The two estimates are consistent, given the 
uncertainties 

Assuming an H^O/H^O ratio of 500, the XCLASS re- 
sult leads to an H 2 column density of 1.5 x 10 18 cm~ 2 . 
With a source diameter of 4" and a distance of 1.7 kpc, 
we derive a total number of water molecules in the dense 
core of 1.22 x 10 52 . Hunter et al. (2006) detected four 
dense cores within NGC 6334 I, of which two are hot 
(T > 75 K). Only these two cores may be emitting in 
the high energy water lines. The derived range of masses 
depends on the assumed dust temperature. Excluding 
scenarios cooler than 50 K, i.e., cooler than the excita- 
tion temperature of H3, 7 and HDO, the total mass of the 
hot cores is between 20 and 60 M Q , which yields a water 
abundance ofl.7xl0~ 6 with a factor of 2 accuracy. This 
water abundance is about two orders of magnitude below 
the value found in some other hot cores (e.g., Herpin et 
al. 2012). The measured abundance ratio of H^O/H^O 
is 3.75 ± 1, in good agreement with the 18 0/ 17 isotopic 
ratio (Wilson & Rood 1994). HDO is a factor of 2 less 
abundant than H3, O, and therefore has abundance simi- 
lar to that in the outflow, a factor of ten higher than the 
cosmic D/H ratio. 

4. RADIATIVE TRANSFER MODELS 

To model H^O lines, we use the spherically symmetric 
Monte Carlo radiative transfer code RATRAN (Hoger- 
heijde & Van der Tak 2000). The source model consists of 
a number of shells with varying H2 density, water abun- 
dance, temperature, turbulent velocity, and expansion 
velocity. From such a source model, RATRAN calcu- 
lates a grid of spectra integrated along a pencil beam, 
which has to be convolved to the required spatial resolu- 
tion. The molecular data (Einstein coefficients and col- 
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TABLE 5 

Line parameters of the H3, 6 lines associated with the dense core. 



Transition 


T B 


VLSR 


Av 


E U p 






(K) 


(kms" 1 ) 


(kms" 1 ) 


(K) 


422 


- 4i 3 


1.9 ±0.4 


-8.0 ±0.1 


5.4 ±0.3 


454.5 


h-2 


-3 3 i 


0.85 ±0.06 


-7.8 ±0.1 


8.1 ±0.2 


454.5 


524 


- 4 3 i 


0.7 ±0.1 


-8.3 ±0.3 


6.1 ±0.4 


599.1 


532 


- 4 4 i 


0.52 ±0.02 


-8.6 ±0.6 


7.8 ±0.1 


732.4 


(124 


-7 17 


0.09 ±0.02 


-8.0 ±0.2 


6.8 ±0.2 


867.7 


634 


°41 


0.43 ±0.27 


-6.9 ±0.5 


2.9 ±0.5 


934.2 



"Tentative detection. 



TABLE 6 

Line parameters of the rare water isotopologues in the dense core. 



Line 




Frequency 


T B 




Ad 






(GHz) 


(K) 


(kms -1 ) 


(kms -1 ) 


H ls O 2 02 


- In 


994.7 


0.6 ± 0.13 


-8.0 ±0.13 


4.7 ±0.26 


H ls O 2xx 


— 2o2 


745.8 


0.47 ±0.05 


-7.95 ±0.08 


6.3 ±0.22 


H 18 3i2 


- 303 


1095.6 


0.6 ± 0.11 


-8.8 ±0.12 


4.8 ±0.30 


H 17 2xx 


— 2o2 


748.5 


0.39 ±0.05 


-8.4 ±0.13 


6.1 ±0.44 


H 17 3i2 


- 303 


1096.4 


0.36 ±0.09 


-8.4 ±0.14 


4.1 ±0.31 


HDO 2 02 


-loi 


919.3 


0.25 ±0.06 


-8.2 ±0.17 


3.9 ±0.40 


HDO 2u 


-Ixo 


1009.9 


0.29 ±0.10 


-8.3 ±0.20 


4.6 ±0.46 


HDO 2n 


— 2()2 


599.9 


0.24 ±0.03 


-8.5 ±0.13 


6.1 ±0.34 
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TABLE 7 

Excitation temperature and column density of H2, 8 0, 
H^ 7 0, and HDO derived with XCLASS. 



Molecule 


1 ex 


Af 




(K) 


(cm -2 ) 


H^O 


43±3 


(3.0 ±0.7) x 10 15 


H^O 


50±4 


(8.0 ±1.0) x 10 14 


HDO 


50±f 


(4.0 ±1.0) x 10 14 



lisional rates calculated by Faure et al. 2007) are taken 
from the LAMDA database (Schoier et al. 2005). The 
source structure is taken from Rolffs et al. (2011), based 
on APEX line and continuum observations, assuming a 
power-law density structure and solving for the tempera- 
ture distribution self-consistently. This spherically sym- 
metric model does not account for the outflow. Since the 
outflow features are dominant in several H^O spectra, 
we had modified the source model accordingly. 

Simply adding shells with a corresponding expansion 
velocity to the existing source model does not reproduce 
the outflow features appropriately. Such an approach 
fails, because, as shown in paper I and §3.3.2, only 17% 
of the continuum is covered by the optically thick out- 
flow. Furthermore, the outflow is collimated almost along 
the line of sight, whereas a modeling approach with an 
expanding shell would reproduce rather an isotropic out- 
flow. To model the collimated outflow, we thus calcu- 
lated two RATRAN models. The first model does not 
include any outflow features. The second one includes 
an isotropic outflow represented by a number of expand- 
ing shells. Before convolving the grid of spectra with the 
beam size of the observations, we combined these two 
models to simulate the desired geometry. 

The RATRAN model without an outflow was calcu- 
lated using the source model from Rolffs et al. (2011), 
with the water abundance as the only parameter. For 
the model with an outflow, we use the same input pa- 
rameters as for the first model, but add a number of 
expanding shells representing the outflow. Temperature, 
density, and water abundance in these additional shells 
are treated as free parameters. The output of both model 
calculations is a grid of pencil-beam spectra. We combine 
these two grids by taking spectra from the model with an 
outflow for points on the grid less than 3.5" from the cen- 
ter, whereas, we use the results from the model without 
an outflow for points further away. The radius of 3.5" 
was chosen to cover approximately 20% of the contin- 
uum, in agreement with previous results. This combined 
grid was convolved with a Gaussian corresponding to the 
respective telescope beam size. 

An alternative, simpler approach is to take the out- 
put of the radiative transfer model without an outflow 
and add the result of a fourth order polynomial fit of the 
outflow based on the observed spectra, also correcting 
for the difference between the model and observed con- 
tinuum. Fits were done using velocities > — 12 kms -1 
(blue lobe) and < +3 kms -1 (red lobe). This method 
gives better constraints on water abundance in the hot 
core and quiescent envelope, because uncertainties due 
to the emission from the outflow are minimized, but the 
model is not self-consistent. 

Assuming a constant water abundance, the model fails 



to reproduce the observed line profiles, and we are forced 
to vary the abundance with temperature. For gas above 
100 K, we derive a water abundance of 6 x 10 -7 , whereas 
below 100 K the abundance drops to 3 x 10 -9 . The H 2 
abundance below 100 K agrees fairly well with theoreti- 
cal models, as well as results of other investigations (e.g., 
Boonman et al. 2003; Chavarrfaet al. 2010). Doty et al. 

(2002) predicts an H 2 abundance ~ 10~ 7 even at dust 
temperatures above 100 K due to the destruction by ions, 
such as C + , H + , , and He + . However, their model 
does not include formation of water on grain surfaces, or 
continuous sublimation of ices infalling from the cold en- 
velope toward the center of the protostar. Other models, 
such as Ceccarelli et al. (1996) and Rodgers & Charnley 

(2003) suggest a water abundance of 10~ 4 in warm gas, 
which is approximately two orders of magnitude higher 
then our result. However, it is important to note that the 
temperature and the amount of hot gas might be over- 
estimated in the RATRAN model due to the assumed 
geometry. The assumed spherical symmetry results in a 
lower escape probability of the cooling photons, than for 
other geometries. This lower escape probability may lead 
to an overestimation of the temperature in the inner part 
of the star-forming region. The water abundance in the 
outflow cannot be determined, since the outflow model 
is mainly defined by features due to the warm and dense 
outflow, for which no independent mass estimate exists. 
However, the total number of water molecules in the out- 
flow is 1.2 x 10 51 , which results in a total outflow mass of 
4.7 x 10~ 2 M , assuming an H 2 abundance of 4.3x 10~ 5 
as found in the broad outflow (paper I). Thus the mass 
is only about 5% of the mass of the broad outflow as 
determined by Leurini et al. (2006), which may explain 
why this outflow feature is not detected in CO. We also 
compared the model results with observations of H3> 8 
and H^O, using an H^O/H^O abundance ratio of 500 
and an H^O/H^O ratio of 1600. 

The best fit results for our model are show in Figures 8 
and 9. The results of models, where the outflow was only 
fitted and not actually modeled (blue lines) show that 
the main component is reproduced quite well. One ex- 
ception is ortho-H20 3i2 — 22i, which is a factor of three 
stronger than predicted by our model. Another exception 
is the para-H20 624 — 717 transition, which is predicted to 
be very week (0.015 K), with the observed line approxi- 
mately seven times stronger (0.11 K at a velocity of -7.7 
kms -1 ). One explanation could be that the 624 — 7n 
line shows maser activity. Alternatively this line could 
be also blended with other lines, especially considering 
the rather large line width (6.5 kms -1 ) and taking into 
account the richness of the spectrum (see Fig. 8). How- 
ever, we could not identify potentially interfering lines in 
the major molecular spectroscopy databases (JPL, Pick- 
ett et al. 19980; CDMS, Muller et al. 200©. 

For lines with very high optical depths, i.e., ortho-H20 
lio — loi or para-H20 In — 0oo, the model calculations 
predict some emission in the line wings, which is not seen 
in the observed spectra. Modeling the outflow is difficult. 
With the exception of the ortho-H20 lio — loi > where we 
cannot find an acceptable model, there is a qualitative 

4 http://spec.jpl.nasa.gov 

5 http://www.astro.uni-koeln.de/cdms 
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agreement between the model and observed line shapes 
(emission, absorption, or a P-Cygni profile). However, 
there are significant quantitative differences between de- 
tailed line shapes and line strengths predicted by the 
model and the observations, which indicate a limitation 
of the one-dimensional approach. The lines of the rarer 
isotopologues are not well reproduced (Figure 10). The 
strengths of the modeled H 2 8 lines are typically only 
60% of the observed ones, whereas modeled H 2 7 lines 
are much too weak. This indicates that the warm ma- 
terial might reside in small, dense clumps and thus the 
optical depths of the water lines may be underestimated 
by the models. This shows that a much more realistic 
three dimensional radiative transfer model is required to 
fit the water lines properly, which in turn requires de- 
tailed interferometric spectral maps of the source. 

5. SUMMARY 

The H3> 6 abundances in the cold envelope (< 100 K) 
and the foreground clouds of 3 x 10~ 9 and 1.5 x 10~ 8 , re- 
spectively, agree fairly well with model predictions for 
cold regions (e.g., Ceccarelli et al. 1996, Doty et al. 
2002), where freeze-out takes place. Chemical models 
also reproduce the water abundance of 4 x 10~ 5 found 
in the broad outflow component, where we expect water 
ices to be sputtered from the dust grains. However, the 
water abundance we derive in the hot core (~ 10~ 6 ) is 
two orders of magnitude lower than the value predicted 
by chemical models for hot, high-density material. The 
hot core H2O abundance in NGC 6334 I derived here is 
in agreement with HIFI studies of hot corinos in low- 
mass star forming regions (e.g., Visser et al. 2012; Kris- 
tensen et al. 2012), but in contrast to some previous 
studies of water in high-mass star forming regions (Mel- 
nick et al. 2010; Chavarria et al. 2010; Herpin et al. 
2012) that found water abundances of order 10 -4 at tem- 
peratures above 100 K. A high water abundance in the 
warm gas is only expected if all water molecules subli- 
mate from the dust grains and stay as H 2 in the gas 
phase. A possible explanation for the low water abun- 
dances in NGC 6334 I are time dependent effects — water 
molecules may not have had enough time to fully des- 
orb from the dust grains, or are quickly destroyed in the 
gas phase. Another possibility is that the dust temper- 
ature is moderate, near the threshold for the water ice 
sublimation, and therefore a significant fraction of water 
molecules may still reside on the dust grains. In addi- 
tion, the presence of small-scale structures, such as disks 
or fragmented dense cores, may lead to overestimating 
the amount of the warm H2 based on dust continuum 
observations — if the H 2 column density is overestimated, 
the resulting H 2 abundance may be underestimated. 
Lines of the rare water isotopes suggest that the inter- 
nal structure of this source may be more complex and, 
therefore, future, more detailed investigations of the H 2 
abundance have to take into account the detailed struc- 
ture of the embedded cores. 

The two detected outflow components in NGC 6334 I 
are physically very different. The broad outflow, which 
has been previously detected in CO (Leurini et al. 2006), 
is of moderate temperature (~ 50 K) and density (n ~ 
10 4 cm' 3 ). The narrow outflow (T ~ 100 K, n ~ 2 x 10 7 
cm -3 ), which, to our knowledge, has not been identified 



in other tracers, is warmer and denser, which confirms 
that water is an excellent tracer of dense outflows (San- 
tangelo et al. 2012). Whether these two components are 
actually two different outflows, or different regions within 
the same outflow cannot be conclusively answered. The 
fact that the central velocity of both outflow component 
is -10 kms -1 points to the latter interpretation. 

The ortho/para ratio in water is 3 ± 0.5 in the fore- 
ground clouds, as well as the outflow. This is consis- 
tent with the scenario, in which water molecules are 
formed with an ortho/para ratio of three, and the re- 
actions which would cause a spin flip, such as those with 
H + and H^, arc of minor importance. The abundances 
of the rarer water isotopologues are consistent with the 
elemental abundance ratios. The H 2 8 0/H 2 7 ratio in 
the dense core is 3.75 ± 1, in agreement with the ratio of 
3.7 found in the broad outflow component (paper I) and, 
within the uncertainties, consistent with the 18 0/ 17 iso- 
topic ratio. 

HDO is detected in the hot core, as well as the broad 
outflow component. The HDO/H 2 ratio in both re- 
gions is ~ 2 x 10~ 4 , a factor ten above the cosmic abun- 
dance of deuterium, and close to the value found in Earth 
oceans. However, this ratio is much lower than that typ- 
ically found in low-mass star forming regions (~ 0.04, 
Codclla et al. 2010), and at the lower end of values 
found in other high-mass star forming regions (Bergin 
et al. 2010; Comito et al. 2003, 2010; Van der Tak et 
al. 2006). HDO is formed in excess in cold clouds, cither 
in the gas-phase or on dust grains, where many deuter- 
ated molecules are abundant, whereas at higher tempera- 
tures (> 50 K) HDO/H 2 ratios similar to the one found 
here are predicted (Roberts & Millar 2000). The high 
HDO abundance in Orion KL was interpreted as mate- 
rial recently desorbed from dust grains, and thus the high 
HDO/H 2 ratio reflects the physical conditions from the 
past. Thus the low HDO/H 2 ratio in NGC 6334 I, 
in agreement with the relatively low water abundance, 
is another indication that most water, especially water 
formed during the cold early stages of star-formation, is 
still frozen-out onto dust grains. 
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Fig. 8. — Comparison of the radiative transfer model results with observations of para-H^O (black). Blue lines represent models where 
the outflow component was fitted by a polynomial function, whereas the red lines show results for calculation where the outflow was 
incorporated into the model. The intensity scale is brightness temperature, in K, corrected for the aperture efficiency. 
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Fig. 10. — Comparison of the radiative transfer model results (red) with observations (black) for HJ^O and H^O. 



van der Wiel, M. H. D., van der Tak, F. F. S., Lis, D. C, et al. 

2010, A&A, 521, 53 
van Dishoeck, E. F. & Helmich, F. P. 1996, A&A, 315, 177 
van Dishoeck, E. F., Kristensen, L. E., Benz, A. O., et al. 2011, 

PASP, 123, 138 



Visser, R., Kristensen, L. E., Bruderer, S., et al. 2012, A&A, 537, 
A35 

Walsh, A. J., Longmore, S. N., Thorwirth, S., et al. 2007, 

MNRAS, 382, 35 
Willner, S. P., et al. 1982, ApJ, 253, 174 
Wilson, T. L. & Rood, R. 1994, ARA&A, 32, 191 



